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Abstract 

We propose to test CP symmetry in the decay J ftp — > 7 </></>, for which large data sample exists 
at BESII, and a data sample of 10 10 J/tfj's will be collected with BESIII and CLEO-C program. 
We suggest some CP asymmetries in this decay mode for CP test. Assuming that CP violation is 
introduced by the electric- and chromo-dipole moment of charm quark, these CP asymmetries can be 
predicted by using valence quark models. Our work shows a possible way to get information about 
the electric- and chromo-dipole moment of charm quark, which is little known. Our results show that 
with the current data sample of J/i/j, electric- and chromo-dipole moment can be probed at order of 
10 _13 e cm. In the near future with a 10 10 data sample, these moments can be probed at order of 
10~ 14 e cm. 

PACS numbers: 11.30.Er,13.25.Gv 
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CP violation was first discovered in if-system in the last century |IJ. At the beginning of this cen- 
tury, CP violation was also found in 5-system thanks to two 2?-factoriesl.2 a . Except these, there is no 
experimental indication for CP violation. With our current knowledge interactions between elementary 
particles are described by the Standard Model(SM). In SM CP symmetry is violated by complex CKM 
matrix elements. Although SM is successful and can explain CP violation effects found in K- and B- 
system, it is generally believed that SM is only an effective theory for energy scale around lOOGeV. It 
is expected that new physics beyond SM will emerge. Therefore, it is important to find CP violation 
somewhere else, where SM predicts no or tiny CP violation, because this can be an indication for new 
physics. 

It is expected that effect of CP violation is very small in general. Hence, Experimental study of CP 
violation needs large data samples. With progress at BEPC and CLEO-C large data samples of J/ip 
decays exist or will be collected. At BEPC with BESII detector there are already 5.8 x 10 7 J/ip-dec&y 
events. With upgraded BEPC 3 and CLEO-C programjlj, a data sample of 10 9 ~ 10 10 J/ip will be 
collected. These huge data samples are very suited for CP test. But, not every decay mode of J/ip can 
be used for CP test. For J/ip decay into a particle and its antiparticle, CP test is not possible if these 
particles are spinless or their polarizations are not observed [JJ]. CP test is only possible if polarizations 
of decay products are measured. The decay J/ip — > AA is an example, where the polarizations can be 
measured though subsequential decays of A and A. CP test with this decay mode has been studied in 
CP test can be made for three-body decays, even without knowing polarizations of decay products. In this 
letter we propose to study CP test in the decay J/ip — » 700. This decay can provide useful information 
about electric dipole moment (EDM) of charm quark. The reason for choosing this channel is because 
is a very narrow resonance just above KK threshold and can be clearly identified by its K + K~ decay 
mode in experiment. In principle J/ip — ► ^ypp could also serve the purpose, but experimentally the broad 
width of p meson makes it impossible to get clean sample for the channel. By using large charmonia 
data samples it has been also discussed in [7] about possible CP test, e.g., ip' — » J/ipirir, and in [Hj about 
parity test in the decay mode J/ip — ► ^£ + £~. 

From SM EDM of quarks and leptons are very small(See Reviews in and references therein). If 

EDM of quark is found to be nonzero, it is likely an indication for new physics. Because the operator for 
EDM does not converse helicities of quarks, its effect is suppressed in a high energy process by a factor 
rriq/E, where m q is the quark mass and E is a large energy scale. For light quarks, useful information can 
be obtained through measurement of EDM of neutron j^j. So far there is no experimental information 
about EDM's of heavy quarks, like charm- and bottom quark. J/ip decays can provide information of 
EDM of charm quark and has the advantage that the effect of EDM will be not suppressed, because the 
large energy scale is around m c . Since in the decay a cc pair is annihilated into a photon and gluons, it 
also provides a way to detecting chromodipole moment of charm quark. In this letter we assume that CP 
violation is introduced by electric- and chromo-dipole moment and study their effect in the decay. The 
moments are defined by effective Lagrangian: 

Lcp = -i^H^F^c - i^c l5 a^G^c, (1) 

where d c is the electric dipole moment, d c is the chromodipole moment. With the decay mode, it is also 
possible to detect the electric dipole moment of s-quark. But its effect is suppressed by m s /m c . In this 
letter we will only consider dipole moments of c-quark. 

Unfortunately, predictions for the decay J/ip — > 700 can not be made without models, especially, 
because there is no reliable method in theory how to determine the formation of bound states like 0. 
For J/ip a good approximation is to use a nonrelativistic wave- function for describing J/ip as a bound 
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state of a cc pair, because c-quark can be taken as a heavy quark. In this work we will also take this 
approximation for (j) meson by taking it as a nonrelativistic bound state of a ss quark pair, where the 
s-quark should be taken as a constituent quark, but not as a current quark. The constituent s-quark mass 
should be taken as m s ~ This approximation was employed in radiative decay into a hadronpT]. 

In this approximation, the decay amplitude is factorized into a partonic part which is for the process 
cc — > 7SSSS, and a nonperturbative part which is related to wave functions of hadrons at the origin. 
The partonic part can be calculated with perturbative theory. The wave functions at the origin can 
be determined by experiment. In this letter we will take this approximation. It should be noted that 
with this approximation reliable predictions for various distributions may not be obtained. But one may 
expect that for our integrated observables defined in below it can be a good approximation. Especially, 
our integrated observables will not depend on wave functions at the origin. 
We consider the decay in the rest-frame of J/ip 

J/i/j(P) - 7 (fc) + 0(p x ) + ^(pa), (2) 

where momenta are given in brackets. Because the two <j> mesons are identical particles, we require 
Pi > p\ to distinguish them in experiment. The decay amplitude can be written as: 

T(J/V ->7#) = ^(Pi»P2)£t, (3) 

where £{(i = 1,2,3) is the spacial component of the polarization vector of J/tp. With the amplitude, a 
density matrix for the decay can be defined as: 

%(P1,P 2 ) = J] A i (pi, P2 )yl*(pi,p 2 ). (4) 

spin 

The sum is over spins of all particles in the final state. This matrix can be decomposed in terms of tensor 
built with pi, P2, dij and £ijk- A similar decomposition can be found in 0. The density matrix contains 
all information about the decay. Any experimental observable can be predicted if we know the density 
matrix for production. We consider the situation where J ftp is produced through e + e~ annihilation in 
their center-mass frame, as at BEPC or CLEO-C. The density matrix pij for the production takes the 
form: 

1 1 r - 1 
pij = -6 ij --(k + ik +j --5 ij ), (5) 

where k + denotes the moving direction of e + . With these density matrices, the expectation value of any 
experimental observable O constructed with k + , pi and p 2 can be predicted by 

(°) = Tfl dT{2irf5(P -V1-V2- QOfHjRji, (6) 

where dT is the differential phase-space of 7^ with p® > p®, M is a normalization factor so that (1) = 1. 
Symmetries will constrain the form of the density matrix. If CP symmetry holds, we have 

Rij (Pi, P2) = Rij (-pi, -P2)- (7) 

For CP test one can check this relation to see if it holds. A convenient way to check this is to use some 
integrated observables. We find two CP odd observables which we are particularly interested in: 

Ox = k+ • pik+ • (p! x p 2 ), 2 = k+ • p 2 k+ • (pT x p 2 ), (8) 
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where pi or p2 are the direction of pi or P2 respectively. From these oberservables, One can define 
CP -asymmetry as 

B i = {9{O i )-9{-O i )) (i = 1,2), (9) 

where 6(x) = 1 if x > and is zero if x < 0. If these asymmetries are not zero, CP symmetry is violated. 
In this letter, we will make predictions for these asymmetries by taking CP violating interactions as given 
in Eq.(l). 
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Figure 1: A typical Feynman diagram for J ftp — > jcpcp, other diagrams are obtained by changing at- 
tachments of gluons to different places. Replacing ccg and eery coupling with the corresponding dipole 
interaction one obtains CP violating amplitude. 



As discussed before, we will use nonrelativistic wave-functions for bound states. A typical Feynman 
diagram is given in Fig.l. From which one can read out the contribution to the scattering amplitude for 

c (li) + c(q 2 ) — > 7 (A;) + s(k{) + s(k 2 ) + s(k^) + s(k<i), which takes the form: 

T 1 (c(q 1 )+c(q 2 ) ^ 7 (fc)+s(fci) + 5(A;2)+s(fc3) + «(fc4)) = «(92)C a6 '^u( g i)ii(fci)7,n(fc4)«(fe)7^ 6 Kfc2), 

(10) 

where C ab) ^ l/ is a matrix with Dirac- and color indices. To obtain the contribution from Fig.l to the 
hadronic decay amplitude, one needs to project the quark pairs into hadron states. The projection is 
made by the replacement 

Ui{qi)vj{q 2 ) -> 3(2 J p ^-W (°) ^7 ' Ql + m c)l ' £(■**) (~7 ' g2 + m c )] fj - , 
v l {k 2 )u j (k 1 ) -> — — ^ 3/2 ^(0) [(~7 • k 2 + m s )7 • e*(pi)(7 • fci + m,)]^ , 
Vi(k4)uj(k3) -> g72^-^^(0)[(-7-fc4 + m s )7-e*(p2)(7-fe + m a )] yj (11) 

where stand for Dirac- and color indices, s(P), s(j>i) and e(p2) is the polarization vector of J/ip(P), 
4>{p\) and 4>(p 2 ), respectively. V\7/i/>(Vv>) is the wave function of J/^{4>). In the projection momenta of 
quarks and masses should be taken as: 

P M J/V , 

91 = 92 = 77- m c = n , 
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It should be noted that from the amplitude in Eq.(lO) there is another contribution to the hadronic 
amplitude which is obtained by exchanging p\ <-» p2- The calculation is straightforward. The total decay 
amplitude is too complicated to be given here. A numerical program for the amplitude can be obtained 
from authors. 

With the approximation explained in detail, numerical results can be obtained. We use a s (m c ) = 0.28 
and masses which are fixed by experimental values of hadrons given in Eq.(12). We obtain the decay 
width: 

r(J/V - 700) = 1-83 x 10- 3 |^(0)| 2 |^(0)| 4 GeV- 8 . (13) 

If we use the value of \ipj/^(0)\ 2 = 0.055GeV 3 , determined from study with potential models^! and 
lattice QCD^I], and experimental results of T to t = 87KeV and Br(J/ip — » 700) = 4.0 x 10 -4 , we can 
obtain 

|VV,(0)| 2 « 0.019GeV 3 . (14) 

It should be noted that in our approximation the CP asymmetries in Eq.(9) do not depend on these wave 
functions at the origin. With the same parameters we obtain: 
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A statistic sensitivity to d c and d c can be determined from these results by requiring that the asymmetry 
generated by these dipole moments should be larger than the statistical error. The statistical error for 
these asymmetries is given by 
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where N events is the number of the decay events. With the 5.8 x 10 7 J/tp data sample at BES, the 
sensitivities of our CP asymmetries to these dipole moments are 



dr 



1.4 x 10~ 13 e cm, d c ~ 4.5 x 10 _13 e 



cm. 



With a 10 10 data sample which will be collected in the near future, the sensitivities are: 



d r 



1.2 x 10" 14 e cm, d c ~ 3.6 x 10~ 14 e 



cm. 



(17) 
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It is instructive to compare these sensitivities to the upper bound of EDM of neutron [To] and 
electron|16| from experiment. The upper bounds read: 



d n < 6.3 x 10 
d e < 4.3 x 10 



-26 



e cm, 



-27 



e cm. 



(19) 



The implication of these upper bounds for physics beyond SM has been studied extensively (See e.g., |17j 
and references therein). Predictions from any reasonable extension of SM for EDM of light quarks and 



electron should not be larger than these upper bounds. Comparing these upper bounds, the sensitivities 
are too large. However, in some extensions of SM predicted EDM of a quark is proportional to the mass 
of the quark(e.g., see |181 I19j ). It is possible to obtain EDM of c-quark which is comparable to these 
sensitivities. Hence, CP test proposed in this letter will provide information for possible new physics 
beyond SM. 

To summarize: In this letter we have proposed to test CP symmetry in the decay J ftp —* jcpcp by 
using large data samples of J/ tp which is already collected at BES and will be collected with BESIII and 
CLEO-C program. We introduce some CP asymmetries for CP test with this decay mode. Assuming 
that CP violation is introduced by the electric- and chromo-dipole moment of charm quark, we give 
predictions for these asymmetries by using wave functions for hadrons. Experimental information is 
available for electric dipole moments of light quarks through measurement of the electric dipole moment 
of neutron, but there is no information for electric dipole moments of heavy quarks. Our work shows a 
possible way to obtain this information. With the current data sample of J/tp at BES, the electric- and 
chromo-dipole moment can be probed at order of 10 -13 e cm. In the near future with a 10 10 data sample, 
these moments can be probed at order of 10 _14 e cm. 
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